The effect of the hydrogen coverage on the Young modulus of graphene 
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We blend together continuum elasticity and first principles calculations to measure by a computer experiment 
the Young modulus of hydrogenated graphene. We provide evidence that hydrogenation generally leads to a 
much smaller longitudinal extension upon loading than in pristine graphene. Furthermore, the Young modulus 
is found to depend upon the loading direction for some specific conformers, characterized by an anisotropic 
linear elastic behavior. 
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PACS numbers: 81.05.ue, 62.25.-g, 71.15.Nc 

The hydrogenated form of graphene ! (also referred to as 
graphane) has been systematically investigated by Wen et al. , 
proving that in fact there exist eight graphane isomers. They 
all correspond to covalently bonded hydrocarbons with a C:H 
ratio of 1 . Interesting enough, four isomers have been found 
to be more stable than benzene, indeed an intriguing issue. By 
hydrogenation, a set of two dimensional materials is generated 
with new physico-chemical properties: for instance, graphane 
is an insulator ' , with an energy gap as large as ~ 6 eV , 
while graphene is a highly conductive semi-metal. Further- 
more, disordered hydrogenated samples show electronic and 
phonon properties quite different than the pristine material . 

As far as the elastic behavior is concerned, it has been 
proved that hydrogenation largely affects the elastic moduli 
as well. By blending together continuum elasticity theory 
and first principles calculations, Cadelano et al. have deter- 
mined the linear and non linear elastic moduli of three sta- 
ble graphane isomers, namely : chair- (C-), boat- (B-), and 
washboard- (W-) graphane. The resulting picture is very in- 
teresting: B-graphene is found to have a small and negative 
Poisson ratio, while, due to the lack of isotropy, C-graphane 
admits both softening and hardening non linear hyperelastic- 
ity, depending on the direction of applied load. 

Although full hydrogen coverage is possible both in ordered 
and disordered configurations, 5 it is more likely that a typical 
experimental processing procedure generates samples with a 
C:H ratio larger than 1 . In other words, we must admit that 
graphane could exist not only in a large variety of conform- 
ers, included the amorphous one (hereafter referred to as A- 
graphane), but also in several forms characterized by different 
stoichiometry. 

In this work we investigate the variation of the Young mod- 
ulus of graphane versus the hydrogen coverage, in three stable 
conformers and in the amorphous one. The work is addressed 
to establish whether an incomplete sp 3 hybridization affects 
the elastic behavior and which is the trend (if any) of varia- 
tion of the Young modulus versus hybridization, an issue of 
large technological impact. Our approach combines contin- 
uum elasticity (used to define the deformation protocol aimed 
at determining the elastic energy density of the investigated 
systems) and first principles atomistic calculations (used to 
actually calculate such an energy density and the correspond- 
ing elastic moduli). 

Atomistic calculations have been performed by Density 



Functional Theory (DFT) as implemented in the QUAN- 
TUM ESPRESSO package . The exchange correlation po- 
tential was evaluated through the generalized gradient ap- 
proximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) 
parameterization , using Rabe-Rappe-Kaxiras-Joannopoulos 
(RRKJ) ultrasoft pseudopotentials ' n . A plane wave basis 
set with kinetic energy cutoff as high as 24 Ry was used 
and the Brillouin zone (BZ) has been sampled by means of 
a (4x4x1) Monkhorst-Pack grid. The atomic positions of the 
investigated samples have been optimized by using damped 
dynamics and periodically-repeated simulation cells. Accord- 
ingly, the interactions between adjacent atomic sheets in the 
supercell geometry were hindered by a large spacing greater 
than 10 A. 

The Young modulus of the structures under consideration 
has been obtained from energy-vs-strain curves, correspond- 
ing to suitable deformations applied to samples with different 
hydrogen coverage, namely: 25%, 50%, 75%, and 100%. The 
corresponding simulation cell contained 8 carbon atoms and 
2, 4, 6, and 8 hydrogen atoms, respectively. For any possible 
coverage (other than the 100% one), several different geome- 
tries (up to 10) have been considered, by randomly distribut- 
ing hydrogen atoms according to different decoration motifs. 
This implies that all data below are obtained through configu- 
rational averages, a technical issue standing for the robustness 
of the present results. 

For any deformation the magnitude of the strain is repre- 
sented by a single parameter Thus, the energy-vs-strain 
curves have been carefully generated by varying the magni- 
tude of £ in steps of 0.001 up to a maximum strain £ max = 
±0.02. This choice warrants that the linear elastic regime 
was carefully explored. All results have been confirmed by 
checking the stability of the estimated elastic moduli over 
several fitting ranges for each sample. The reliability of the 
present computational set up is proved by the estimated val- 
ues for the Young modulus (E) and the Poisson ratio (v) of 
pristine graphene (corresponding to 0% of hydrogen cover- 
age), respectively 349 Nm -1 and 0.15, which are in excellent 
agreement with recent literature. 7 ' 1 ' -15 Similarly, our results 
E = 219 Nm _1 and v = 0.21 for stoichiometric C-graphane 
(corresponding to 100% hydrogen coverage) agree with data 
reported in Ref. 7, 16, and 17. 

Among the systems here investigated, stoichiometric C- 
graphane is elastically isotropic due to its crystallographic 
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symmetry, while its non stoichiometric conformers with 25%, 
50% and 75% hydrogen coverage (see Fig. 1) are so by 
assumption (which is indeed reasonable by only assuming 
that the hydrogen decoration in real samples is totally ran- 
dom). Similar properties hold for A-graphane. On the other 
hand, stoichiometric B- and W-graphane show an orthorhom- 
bic symmetry, which causes an anisotropic linear elastic be- 
havior. Their non stoichiometric conformers could became 
isotropic due to the achieved disorder degree, as discussed be- 
low. 

The two-dimensional elastic energy density (per unit of 
area) for systems with orthorhombic symmetry (namely, B- 
and W-graphane) is expressed as ' 

U = 4 + + ^12£xv£vv + 2^444 ( 1 ) 

where x (y) indicates the zigzag (armchair) direction in the 
honeycomb carbon lattice and e, 7 (; , j ~xor y) is the infinites- 
imal strain tensor. By applying a loading at an angle 9 with 
respect to the zigzag direction, the corresponding Young mod- 
ulus E(9) is written as 7 

E{B) = ' V 5 r (2) 

W + ^22C 4 + ( * u %* n ~ 2V n ) 

where c = cos 9, and s = sin0. Eqs.(l) and (2) are also valid 
for isotropy systems (namely, C- and A-graphane) provided 
that it is set: = ^22 and 2^44 = — ^12 (Cauchy rela- 
tion). 

The Young modulus E(9) can be directly obtained from 
the linear elastic constants %j, in turn computed through 
energy- vs-strain curves corresponding to suitable homoge- 
neous in-plane deformations. Only two in-plane deformations 
must be applied to obtain all the independent elastic con- 
stants, namely: (i) an uniaxial deformation along the zigzag 
(or armchair) direction; and (ii) an hydrostatic planar defor- 
mation. Nevertheless, for the validation of the isotropicity 
condition, two more in-plane deformations must be further 
applied: (iii) an axial deformation along the armchair (or 
zigzag) direction; and (iv) a shear deformation. The strain 
tensors corresponding to applied deformations depend on the 
unique scalar strain parameter £ 7 ' 15 , so that the elastic en- 
ergy of strained structures defined in Eq.( 1) can be written as 
U(Q = U + 5f/ (2) C 2 + 0(C 3 ), where U is the energy of the 
unstrained configuration. Since the expansion coefficient £/' 2 ' 
is related to the elastic constants %j, a straightforward fit pro- 
vides the full set of linear moduli for all structures. Repeating 
this procedure for several randomly distributed hydrogens in 
different configurational symmetries (corresponding to C-, B-, 
W-, and A-graphane), the complete set of elastic constants %j 
has been evaluated (as an average value) at various coverages. 

The synopsis of the calculated elastic constants %j for all 
graphane samples here investigated is reported in Table I, from 
which quite a few information can be extracted. As a pre- 
liminary remark, we stress that each hydrogenated conformer 
is characterized by a specific hydrogen arrangement and by 
a different buckling of the carbon sublattice. Moreover, in 
non stoichiometric samples hydrogen atoms can migrate (as 
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Figure 1. Pictorial representations of some hydrogen motifs corre- 
sponding to 25% (Panel a), 50% (Panel b), and 75% (Panel c) cov- 
erage. Hydrogen atoms are indicated by red (dark) circles, while hy- 
drogen vacancies by gray (light) circles which are randomly placed 
on the top or bottom of the graphene sheet (black sticks). Shaded 
areas represent the simulation cell. Up to 10 different decorations for 
each conformer have been used in the calculations. 



indeed observed during relaxation) due to the presence of 
nearby unsaturated bonds, thus affecting the resulting symme- 
try of the conformer. These features add further details to an 
already complex situation, introducing another source of dis- 
order. We, therefore, argue that the amorphous configuration 
should be considered the most probable phase. Alternatively, 
the above mentioned migration of hydrogen atoms could in- 
deed generate phase separation into C-H and pure C domains, 
as discussed by Lin et air . 

By inserting the values reported in Table I into the Eq.(2), 
we obtain the Young modulus E(9) for any hydrogen cover- 
age, as shown in Fig. 2. As a general feature, we state that 
the change in hybridization has largely reduced the property 
of longitudinal resistance upon extension, as described by the 
greatly reduced value of the Young modulus with respect to 
ideal graphene (bold dashed line in Fig. 2). Simulations pro- 
vide evidence that this is mainly due to the fact that sp 3 hy- 



Table I. Independent elastic constants (units of Nm~ ) are shown for 
different values of the hydrogen coverage. The 0% coverage corre- 
sponds to the pristine graphene, where ^n=357 ±7, and ^12 52±1 1. 
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Figure 2. Polar plot of the Young modulus E(6) (see Eq.(2)) shown 
as function of the hydrogen coverage. In each quadrant a given con- 
former of graphane is represented, as indicated by labels. The Young 
modulus of pristine graphene is also shown by the bold dashed line. 



bridization creates tetrahedral angles (involving 4 carbons and 
1 hydrogen) which are easily distorted upon loading. In other 
words, softer tetrahedral deformations are observed, rather 
than bond stretching ones as in ideal graphene. Equivalently, 
we can conclude that deformations upon loading are basically 
accommodated by variations of the tetrahedral angles. 

The full picture about the elastic response is provided in 



Fig. 2, by considering the 9 -dependence of the calculated 
Young modulus for each graphane conformer. The corre- 
sponding value of ideal graphene is reported as well for com- 
parison. In such a polar plot, an isotropic behavior corre- 
sponds to a circular line. This is always found for C- and 
A-graphane (as well as for ideal graphene), as expected. On 
the other hand, above 75% hydrogen coverage, the B- and W- 
conformers show a strong 9 -dependence in their elastic be- 
havior (or, equivalently, of their Young modulus). Interesting 
enough, well below this coverage they almost behave isotrop- 
ically. 

Deviations from the isotropic elastic behavior are quantita- 
tively predicted by the calculation of the = 4^44 / \ + 
^22 — 2^12) ratio, which should be 1 for ideally isotropic 
systems. Data reported in Table I provide si = 1.0 ±0.05 
for any C- and the A-conformers, thus confirming their elas- 
tic isotropicity. On the other hand, non stoichiometric B- 
and W-samples with hydrogen coverage smaller than 75% 
display a comparatively small anisotropy, corresponding to 
£/ = 1.0 ±0.2. Finally, these same systems with higher 
hydrogen coverage recover a full anisotropic behavior is 
largely different than 1). 

In conclusion, we have presented and discussed first princi- 
ples calculations showing that the elastic behavior of graphene 
is largely affected by hydrogen absorption by the actual cover- 
age. In particular, the Young modulus is greatly reduced upon 
hydrogenation, as also previously discussed in Ref. 17. An 
incomplete coverage generates a large configurational disor- 
der in the hydrogen sublattice, leading to a larger corrugation 
with respect to highly-symmetric C-graphane. Indeed, such a 
corrugation of the carbon sublattice is a key feature affecting 
the overall elastic behavior. 
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